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ABSTRACT 

In this paper, we present a modijied rariational iteration method for solving Fredholm integro-differential 
equations. This study prorides an analytical approximation to determine the beharior of the solution. Moreover, it proves 
the existence and uniqueness in results and convergence of the solution. Finally, some examples are included to 
demonstrate the validity and applicability of the proposed technique. 
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1. INTRODUCTION 

In recent years, there has been a growing interest in the integro-differential equations, which are a 
combination of differential and integral equations. The nonlinear Fredholm integro-differential equations play an 
important role in many branches of nonlinear functional analysis and their applications in the theory of 
engineering, mechanics, physics, electrostatics, biology, chemistry and economics [1-6,22-30]. In this paper, we 
consider the Fredholm integro-differential equations of the type: 

k i 

£ l (x)Z (i \x) = f(x) + yj ' K(x, t)G(Z(t))dt 

( 1 . 1 ) 

with the initial conditions 

Z (r) (a) =b, r= 0,1,2 ,....,(k- l), (1.2) 

where Z <n (x) is the j derivative of the unknown function Z(x) that will be determined, K(x,t) is the kernel 
of the equation, f(x) and Cj(x) are an analytic function, G is nonlinear function of Z and a, b, y, and b are real finite 
constants. Recently, many authors focus on the development of numerical and analytical techniques for integro- 
differential equations. For instance, we can remember the following works: Abbasbandy and Elyas [2] studied 
some applications on variational iteration method for solving system of nonlinear Volterra integro-differential 
equations, Hamoud and Ghadle [6] applied the hybrid methods for solving nonlinear Volterra-Fredholm integro- 
differential equations, Alao et al. [4] used Adomian decomposition and variational iteration methods for solving 
integro-differential equations, Yang and Hou [21] applied the Laplace decomposition method to solve the 
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fractional integro-differential equations, Mittal and Nigam [16] applied the Adomian decomposition method to 
approximate solutions for lractional integro-differential equations, and Behzadi et al. [5] solved some class of nonlinear 
Volterra-Fredholm integro-differential equations by homotopy analysis method. Moreover, several authors have applied 
the Adomian decomposition method and the variational iteration method to find the approximate solutions of various types 
of integro-differential equations [6, 9-11, 16, 21]. 

The main objective of the present paper is to study the behavior of the solution that can be formally determined by 
semi-analytical approximated method as the variational iteration method. Moreover, we proved the existence, uniqueness 
results and convergence of the solutions of the Fredholm integro-differential equation (1.1). 

2. DESCRIPTION OF THE METHOD 

Some powerful methods have been focusing on the development of more advanced and efficient methods for 
integro-differential equations such as the variational iteration method [2, 4, 13, 15]. We will describe this method in this 
section: 

2.1. Modified Variational Iteration Method (MVIM) 

MVIM is applied to solve a large class of linear and nonlinear problems with approximations converging rapidly 
to exact solutions. The main idea of this method is to construct a correction functional form using general Lagrange 
multipliers. These multipliers should be chosen such that its correction solution is superior to its initial approximation, 
called trial function. It is the best within the flexibility of trial functions. Accordingly, Lagrange multipliers can be 
identified by the variational theory [17-19]. A complete review of He’s variational iteration method is available in [13, 15]. 
The initial approximation can be freely chosen with possible unknowns, which can be determined by imposing 
boundary/initial conditions. To illustrate, we consider the following general differential equation: 

LZ(x)+NZ(x)=f(x), (2.1) 

Where, L is a linear operator, N is a nonlinear operator and/(xj is inhomogeneous term. According to variational 
iteration method [4], the terms of a sequence Z„ are constructed such that this sequence converges to the exact solution. The 
terms Z„ are calculated by a correction functional as follows: 

z„ +1 (x) = Z (x) + £ A(T)(LZ n (t) + Ny(T) - f(r))dT. (2.2) 

The successive approximation Z n (x), n > 0 of the solution Z(x) will be readily obtained upon using the obtained 
Lagrange multiplier, and by using any selective function Zo. The zeroth approximation Z 0 may be selected using any 
function that just satisfies at least the initial and boundary conditions, with X determined, several approximations Z n (x), n > 
0 follow immediately. 

The VIM has been shown to solve effectively, easily and accurately a large class of nonlinear problems with 
approximations converging rapidly to accurate solutions. 

To obtain the approximation solution of IVB (1.1) - (1.2), according to the VIM, the iteration formula (2.2) can 
be written as follows. 

Z n+ M)=Z n (x)+L- l [A(x)[fy j (x)Z"Xx)-f(x)-r\ b K(x,t)G(Z n (t))dt\\, 

n J a 
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where L _1 is the multiple integration operator given as follows: 

/, (•) = ! j" ...J {-)dxdx...dx{k — times). 

To lind the optimal /ix), we proceed as tbllows: 

k i 

SZ n Jx) = SZJx)+ f(x)- K(x,t)G{Z n {t))dt\\ 

j =0 

= SZ n (x) + A(x)SZ (x)- L'[SZ (x)A (,)]. ^2 

From Equation (2.3), the stationary conditions can be obtained as follows, 

A (x) = 0 ,and 1 +A.(jc)l jMf = 0. 

As a result, the Lagrange multipliers can be identified as !(,) = -1 and by substituting in Equation (2.3), the 
tbllowing iteration 1'ormula is obtained. 


ir/(- r )l 'st' (x-a) r 


Z 0 (,) = L'[^] + X 


b , 


, =0 r\ 

k ^ 

z„jx) = z/b-f I E^wz:;’w-/(b-r(L(.v,()G(Z(fp].n > o. 


(2.4) 


The term V '— ’ y is obtained from the initial conditions, g t (x) + 0 . Relation (2.4) will enable us to 
Lo r\ 

determine the components Z (x) recursively for n > 0. Consequently, the approximation solution may be obtained by 
using 

Z(t) = limZ(f). 


3. MAIN RESULTS 


In this section, we shall give an existence and uniqueness results of Equation (1.1), with the initial condition (1.2) 
and prove it [6-8]. We can write the equation (1.1) in the form of: 


Z(x) 


L -'[ISfl] + y 




h + 


Y l [{ 


1 K(x,t)G(ZJt))dt]-!+'[£ Z ul (i)]. 


?Jx) 


ZM) 


(3.1) 


Such that, 


Z7‘ [J ~^—K(x, t)G(Z n (t))dt ] = 

“ ?t( x ) 


rb (x — t) 

[ -- —K(xJ)G(Z (t))dt 

] °k\lg k (x) 


g L -[^] z c. 

£ &(*) 


■w=Zj. 


jx-vT%(t)_ zU) 
k-Vj k (t) 


-Z 1 (t)dt. 


We set. 
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f(x )~| ^(x-a)' 


■b . 


y(x) = L-' [y—] + X 
£(*) 


Before starting and proving the main results, we introduce the following hypotheses: 

(Hl)There exist two constants a and y > 0, j = 0 ,...,k such that, for any Z ,Z e C(J ,□ ). 


|G(Z,))-G(Z 2 ))|<a|z l -Z 2 | 

and 

\d j (Z,) - D j ' (Z 2 )| < |z, - z 2 1 , 

We suppose that the nonlinear terms G(Z(x)) and d'(Z) = (—)Z(*) = ^z , (D' is a derivative operator), 
j — 0,1,..., k , are Lipschitz continuous. 

(H2)We suppose that for all a <t < X < b, and y = 0,1,..., k, : 


y(x- t) k K(x,t) 


y(x — t) k K(x,t) 

k ! £ k (x) 

- P’ 

k\ 


(x-O t- '^.(t) 

(*-l)!£(r) 


(.r-t) t- '^.(f) 

(*-l)! 


<e„ 


(H3)There exist three functions d’ 0 , and y eC(D,D*), the set of all positive function continuous on 
Z) = {(x,f)eDxD :0<f<x<l} such that: 

0, = max \8 , 6' = max D , and r* = max \y \. 

(H^J^P^^) is bounded function for all X in J =[a,b]. 

Theorem 3.1. Assume that (H1)-(H4) hold. If 
0 < y/ = (ad t + ky''6])(b - a) < 1. 

Then, there exists a unique solution Z(x ) e C(J) to IVB (1.1) - (1.2). 

Proof. 

Let and Z, be two different solutions of IVB (1.1) - (1.2), then 


I z - z, 1=1 


■ y(x- 1 )‘ K(x,t) r (x-t)"-'£(t) . 

-[G(Z,) - G(Z )))dt - X -— [D (Z ) - D‘ (Z ))]dt I 

£(*)*! ’ % k (t)(k- 1 )! 




r(.r -0 JT(JC ,0 


(x)k ! 

<(ae t + ky 9, )(b - a) I Z_ - ZJ, 


G(Z)-G(Z ))lrff-£[ 


.(r-f)‘‘s,(0 


y- £(f)(*-l)! 


II D (Z )- D (Z, ))|* 
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we get (1 - \fr) I Z t - Z, l< 0. Since 0 < \j/ < 1, so \Z t -Z 2 \=0. therefore, z, = Z, and the proof is completed. 

Theorem 3.2. If prohlem (1.1) - (1.2) has a unique solution, then the solution Z n (x) obtained from the recursive 
relation (2.4) using VIM converges when 

0 <</) = (ad 2 i ky 6 )(b - a) < 1. 

Proof. We have from equation (2.4): 

Z n Jx) - Z(x) = z n (x) - Z(x) - (/. [ X = '' > 1 / 1,1 - zli ’(*)]] - K(x,t)[G(Z (t)) - G(Z(t))]dt\). 


If we set, (x) = 1, and W +1 (x) = Z n+1 (x) - Z(x), W n (x) = Z n (x) -Z(x) since W n (a) = 0 , then 


W 


J(x) = W n (x) + \' 


b yK(x,t)(x — t) 


k\ 


[G(Z n (t))-G(Z(t))\dt- u 


>■Z^.(t)(x-t) k ~' 


(k-iy. 


[d 1 (Z (0) - D 1 (Z(t))\dt - (W (X) - W n (a)). 


Therefore, 


r*l yK(x,t)(x — t) L | 44 f»i Y%J)(x - t) k 1 | 

IW ,(x)l< I--llW \adt+Y I 1 -lmaxl y. IIW I dt 

k\ ” t (k- 1)! J " 


<1W I f j* adjt + u 9\ I max I y s I 

/=o 

<1W I (ot6» 2 + ky’9\ )(b - a) =\ W I </>. 
Hence, 

IIW +1 11= max v+eJ I W +1 (x) I 
< <f> max Vte/ IW (x) I 
= ^ IIW II. 


Since 0 < (j) < 1, then II W ll->0. So, the series converges and the proof is complete. 

4. NUMERICAL RESULTS 

In this section, we present the semi-analytical technique based on the VIM to solve Fredholm integro-differential 
equations. To show the efficiency of the present method for our problems in comparison with the exact solution, we report 
absolute error. 

Example 1 . Consider the following Fredholm integro-differential equation. 

Z (x) = e -x + xe' + f xZ(s)ds, 

Jo 

with the initial condition 
Z(0) = 0, 

and the exact solution is Z(x) = xe'. 
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Table 1: Numerical Results of the Example 1 


X 

Exact Solution 

MVIM 

Er MV iM 

0.1 

0.1105170918 

0.1105170888 

3.00xl0 y 

0.2 

0.2442805516 

0.2442805506 

1.00xl0 -y 

0.3 

0.4049576424 

0.4049576457 

2.80xl0 -y 

0.4 

0.5967298792 

0.5967298865 

7.30xl0 -y 

0.5 

0.8243606355 

0.8243606264 

9.10xl0 -y 

0.6 

1.0932712800 

1.0932712700 

1.00xl0 -!i 

0.7 

1.4096268950 

1.4096267890 

1.06xl0 -/ 

0.8 

1.7804327420 

1.7804325720 

1.70xl0 -/ 

0.9 

2.2136428000 

2.2136429590 

1.59xl0 -/ 

1.0 

2.7182818280 

2.7182818230 

5.00xl0 -y 


Example 2. Consider the following Fredholm integro-differential equation. 


Z'(x) = 


5 

4 



— s)Z~(s)ds, 


with the initial condition 


Z(0) = 0, 


and the exact solution is 

Z(x) = x. 


Table 2: Numerical Results of the Example 2 


X 

Exact Solution 

MVIM 

Er MV iM 

0.1 

0.10000 

0.09987 

1.3xl0 -4 

0.2 

0.20000 

0.19974 

2.6xl0 -4 

0.3 

0.30000 

0.29962 

3.8xl0 -4 

0.4 

0.40000 

0.39951 

4.9xl0 -4 

0.5 

0.50000 

0.49941 

5.9xl0 -4 

0.6 

0.60000 

0.59934 

6.6xl0 -4 

0.7 

0.70000 

0.69928 

7.2xl0 -4 

0.8 

0.80000 

0.79925 

7.5xl0 -4 

0.9 

0.90000 

0.89926 

7.4xl0 -4 

1.0 

1.00000 

0.99929 

7.1xl0 -4 


5. CONCLUSIONS 

We present the MVIM for solving Fredholm integro-differential equations. From the computational viewpoint, 
the MVIM is more efficient, convenient and easy to use. The method is very powerful and efficient in finding analytical as 
well as numerical solutions for wide classes of linear and nonlinear Fredholm integro-differential equations. Moreover, we 
proved the existence and uniqueness results and convergence of the solution. The convergence theorem and the numerical 
results establish the precision and efficiency of the proposed technique. 
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